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Abstract 



We study next-to-leading corrections to the integral kernel of the BFKL equation for 
high energy cross-sections in QCD and in supersymmetric gauge theories. The eigenvalue 
of the BFKL kernel is calculated in an analytic form as a function of the anomalous di- 
mension 7 of the local gauge-invariant operators and their conformal spin n. For the case 
of an extended N = 4 SUSY the kernel is significantly simplified. In particular, the terms 
non-analytic in n are canceled. We discuss the relation between the DGLAP and BFKL 
equations in the N = 4 model. 
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1 Introduction 



The Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation |], |2|] is used now together with 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation || for a theoretical 
description of structure functions of the deep-inelastic ep scattering at small values of the 
Bjorken variable x. In this kinematical region the structure functions are measured by 
the HI and ZEUS collaborations || at DESY. For the DGLAP equation the radiative 
corrections to the splitting kernels are well known || . Although the BFKL equation in the 
leading logarithmic approximation (LLA) was constructed many years ago, the calculation 
of the next-to-leading corrections to its kernel was started only in 1989 [H] and completed 
recently [|7], §]. The total cross-section for the colourless particle scattering is expressed 
in terms of the solution G^ilf , q ) of the BFKL equation for the t-channel partial wave 
fj(~(f~c?;t) (J = 1 + w, t = 0) describing the scattering of the reggeized gluons having 

their initial and final transverse momenta if and ~q correspondingly. 

For unpolarized colourless particle collisions one can average the BFKL kernel over 
the angle between the two-dimensional vectors if, q . The eigenvalue of the averaged 
kernel was calculated in 0. The high energy behaviour of the total cross-sections a ~ 
s A is related with the pomeron intercept A. It turns out, that in MS'-renormalization 
scheme the NLO corrections to A are negative and big. However, with the use of non- 
abelian renormalization schemes and the BLM procedure for the optimal scale setting 
it was shown ||, that the BFKL predictions for the total cross-section of the process 
7*7* — > hadrons turn out to be in an agreement with the experimental data obtained by 
the L3 collaboration [ID|. There are other resummation methods (see |TT[ and references 



therein) which lead to quite similar numerical results. 

Because the initial virtual photons of the process 7*7* — > hadrons can be polarized, 
one should know also the NLO corrections to the spin correlators at high energies. It 
is one of the reasons, why in this paper we calculate next-to-leading corrections to the 
intercept A of the BFKL pomeron with an arbitrary conformal spin n. The case n = 
corresponds to the cross-section cr averaged over the angle "d between the photon spins 
(see J7[) and the case n = 2 describes the high energy behaviour s A2 of the contribution 
<72 proportional to cos 2 $ — |. 

On the other hand, the contribution of higher conformal spins is interesting from a 
pure theoretical point of view, because it is related to remarkable mathematical properties 
of the BFKL equation. To begin with, this equation in LLA is invariant under the Mobius 
transformation of holomorphic impact parameters 

apk + b 
cp k + d 

for arbitrary complex values of a, b, c and d [Q. This invariance allows us to find its 
solution for arbitrary momentum transfers q = providing that we know the eigen- 
values A(m, m) of the kernel for all conformal spins n. The conformal weights m and rh 
are related to the anomalous dimension 7 = | + iv of the local operators belonging to a 
basic series of unitary representations of the Mobius group as follows 
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Tl Tl 

m = 7 + -,m = 7--. (2) 

Further, in LLA the homogeneous BFKL equation in the impact parameter representation 
can be written as the Schrodinger equation 

q 2 

H 12 m = E 12 ^> , A = min E 12 , (3) 

where A is the BFKL intercept. The pair Hamiltonian H k i for the interaction of two 
reggeized gluons with their holomorphic coordinates p k , pi and momenta p k = id k , pi = idi 
can be presented as a sum of the holomorphic and anti-holomorphic Hamiltonians JT5[: 



H H = h kl + h* kl . (4) 

Here 

h k i = ln(pfeP/) + — \n(p M )p k + — \n(p k i)pi + 27 , (5) 
Pk Pi 

7 is the Euler constant and p^i = Pk — Pi- 

The holomorphic separability of the Hamiltonian is very important for the solution of 
the generalized BFKL equation [14| for a composite state of n reggeized gluons: 

= A n = -^E. (6) 



The Hamiltonian H in the multi-colour QCD N c — > 00 can be presented as follows |T3 



H=l{h + h*), h = J2 h k)k+1 , [h, h*)=0. (7) 
z k=i 

Due to the commutativity of h and h* the function \l/ has the property of the holomorphic 
factorization |i~3 | 

V = Y1 CrsA^s, (8) 
r<s 

where ip r and if) s are solutions of the Schrodinger equations 

hip r = eip r , h*4> s = ei> s , E = -(e + e) (9) 

with the same values of energies e and e in the holomorphic and anti-holomorphic sub- 
spaces correspondingly. The coefficients C rs should be chosen from the requirement of a 
single-valuedness of ^ in the two-dimensional impact parameter space. 

The holomorphic Hamiltonian h has the integrals of motion q r (r = l,2...n) 



Qr y t Piii2Pi2i3' "Piril PilPi2" -Pir (^0) 



11<12<...<V 

with the properties 



[q r ,q s }=0, [q r ,h] = 0. (11) 
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Moreover, it coincides with the local Hamiltonian of an integrable Heisenberg spin model 



0. 



The integrability of the model is related to a duality symmetry ]16| of h and q r under 
the canonical transformation 

Pr -> Pr,r+1 ~> Pr+1 (12) 

combined with the transposition of these operators. 

In particular, in the case of the compound state of three reggeized gluons, where we 
can use the conformal anzatz 



* ~ = ( -p^) n ( 

^ \P20P30J UoPW 



/ *N P12P30 / 1Q X 

<Pm,m( X , X ) , X = , ( 13 ) 

P10P32 



the duality equation for (p m ~(x, x*) can be written in the pseudo-differential form 

Qm,mfm,m = a ™ a m Vm,m = l^ m l ^l-m,l-m > °"m = — x ) (^) + ■ (^) 

Using this equation, new three-gluon states having the pomeron and odderon quantum 
numbers and the intercept A = were constructed ||17|. Note, that the effective action 



describing the interaction of reggeized gluons with quarks and gluons was derived in ref. 







To investigate a possibility of the holomorphic separability of the BFKL kernel at the 
next-to-leading approximation in QCD and in extended supersymmetric gauge theories 
is another reason for the calculation of the n-dependence of the pomeron intercept. A 
presumable integrability of the BFKL dynamics in the N = 4 supersymmetric field the- 
ory at N c — > oo could be related to remarkable properties of the renormalization group 
equations in this model. It is well known, that the /3-function here is zero in all orders of 
the perturbation theory. Moreover, if we consider conformally invariant evolution equa- 



tions for matrix elements of the quasi-partonic operators introduced in ref. it turns 
out [p0| , that for N = 4 the pair anomalous dimensions for the supersymmetric twist-2 
partonic operators with the Lorentz spin j are proportional to — 1) — It means, 
that for N c — > oo the corresponding Hamiltonian coincides with the local Hamiltonian 
for an integrable Heisenberg spin model similar to that obtained for the reggeized gluons 
(see [^, 15|). The integrability of the evolution equations at the N = 4 theory in the 



multi-colour limit is important for a verification of the Maldacena guess |2T[ 



In this paper we extend the analysis of article to the NLO level. Namely, we calcu- 
late eigenvalues of the angle-dependent eigenfunctions of the NLO BFKL equation for the 
t-channel partial waves in the framework of QCD and in gauge models with an extended 
super symmetry. In the case of QCD this analysis allows us to find main NLO corrections 
to an angle-depended part of cross-sections. On the basis of the study of gauge models 
with an extended supersymmetry one can attempt to find a possible source for the viola- 
tion of the Mobius symmetry, which is important for the calculation of the t-dependent 
BFKL-kernel in QCD at the NLO level. 

The article is organized as follows. In Section 2 we review the results of ref. 0, 
calculate the eigenvalues of the NLO BFKL equation in QCD and discuss the asymptotic 
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behaviour of cross-sections at large s. Section 3 contains the analysis of the gluino and 
scalar particle contributions and our results for eigenvalues of the NLO BFKL equation in 
supersymmetric theories. A summary is given in Conclusion. In Appendix A we discuss 
the reasons for the appearance of the special function Ls3(x) in our formulae and in 
Appendix B we consider the Pomeranchuck singularity in QED and in its supersymmetric 
extension. 



2 NLO corrections to BFKL kernel in QCD 

To begin with, we review shortly the results of ref . |7|] , where the QCD radiative correc- 
tions to the BFKL integral kernel at t — were calculated. We discuss only the formulae 
important for our analysis. 

The total cross-section a(s) for the high energy scattering of colourless particles A, B 
written in terms of their impact factors <&i(qi) and the i-channel partial wave G w (q, q') for 
the gluon-gluon scattering is 



a(s) 



d 2 qd 2 q' 
(2tt)V q' 2 



$A(q) ^B(q') 



a + io ° du 

-ioo llXl 



G u 



So 



\Q\\Q 



(15) 



Here q and q' are transverse momenta^ of virtual gluons and s = 2paPb is the squared 

invariant mass for the colliding particle momenta pa and pa- 
using the dimensional regularization in the MS'-scheme to remove ultraviolet and 

infrared divergences in intermediate expressions, the BFKL equation for G w (q,q') can be 

written in the following form 



uG u ,{q,q 1 ) = 5 D - 2 (q- qi )+ I d D ~ 2 q 2 K(q, q 2 ) G^(q 2 , q x ) 



iD-2, 



(16) 



where 



#(91,92) = 2w(gi)5 D 2 (q x - q 2 ) + K r (q 1 ,q 2 



(17) 



and the space-time dimension D = 4 + 2e for e — > 0. The gluon Regge trajectory u(q) and 
the integral kernel K r (qi,q 2 ) related to the real particle production have been calculated 
in 0, |22|, |2|, |24| and |25|. The final results (see §) can be written in the following 
form correspondingly for the Regge trajectory 




(18) 



3 To simplify equations hereafter (except Appendix B) we omit arrows in the notation of transverse 
momenta if, q , q{. q2, ... , i.e. in our formulae the momenta if, q , qt, gj, ... will be represented 
as g, q' ', gi, q2, ... , respectively. Note, however, that the momenta pa and pb are £>-space momenta. 
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and for the kernel 



K r (qi,q 2 ) 



7^7 ,,-2e 



+a 



+ 



11 2n f \ 1 
Y ~ 3N~J e 

(gi-g 2 ) 2 \V67 



7ri+*r(l-e) i( gi -g 2 ) 2 



In 



9 S V ; 9£ 



1 - C(2)^ 

404 , J/t/n v 56 n/ 
17 + 14 « 3) + 27^ 



jy-3 J ^a„2„2 \ „2 „2 I „2 „2 J I „2 



16gfg| 



Qi Q2 \(l2 Qi 



3+1 + 



i 



(gf + g 2 2 ) 2 _ 2g 2 g 2 2 - 3gf - 
%q\q\ 



+ 



(<? 2 - ? 2 ] 



(qi - q2) 2 (qi + q2) 2 

-d 



In 

(<?? - 4 



16qfq% 
q i\\ n ( ^iliili ~ 12 



(Q1Q2) 



00 dx 



q\ + xg| 



?2, 



(g? + <? 2 2 ) 4 

1 /-CO 



-2Li 2 I -g I +2Li 



In 



1 + x 



1 — X 



(ix j / xq\ 

(?i - <?2) 2 (<?i + <?2) 2 y Wo A y (g 2 - a;gi) 2 V g 2 , 



where Li 2 (x) and ((n) are the Euler dilogarithm and Riemann ^-functions: 

r x dt 

k=l 



Here we introduced also 



r x at 

Li 2 (x) = - / -ln(l-t), CH = E fc " n - 
_ _ s 2 iV c r(l - e) 



(19) 



(20) 



(21) 



" (4tt) 2 + £ 

for the colour group SU(N C ) and gr^ is the QCD coupling constant fixed at the normal- 
ization point /i in the MS'-scheme. 

As it was shown in 0, a complete and orthogonal set of eigenfunctions of the homo- 
geneous BFKL equation in LLA is 



2 \ 7-1 



And 



[22) 



The BFKL kernel in this representation is diagonalized up to the effects related with 
the running coupling constant a s (q 2 ): 



a s {q 2 )N c 



71 



a s {q 2 )N c 

47T 



(23) 



The calculation of the correction 8{n, 7) for all conformal spins is one of main results of 
our paper. 
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2.1 The technique of calculations 



The integration over the angle 6 of the transverse gluon momentum has the form 
r d6e me = 2 dOcosinO) = 2 (* d0T n (cos 9) , 

J-TT JO JO 

where T n (cos#) are the Chebyshev polynomials (below n is a positive integer number). 

In the framework of the dimensional regularization a natural extension of the Cheby- 
shev polynomial is the Gegenbauer polynomial C^(cos#) (see ||26||-||29||) with the index A 
related to the transverse space dimension D — 2 as A = (D — 2)/2 — 1 = D/2 — 2. When 
A = we return to the Chebyshev polynomials because 



C°(cos#) = lim-C*(cos0) = -T n (cos0). 
a^o A n 

To calculate the terms singular at e — > it is convenient to to use the formulae from 
27], E3 for the traceless products 



I n / 2 ] riU — 1 W(n — r> 4- W 
The traceless products are related to the Gegenbauer polynomials C*(cos#) as follows 

{qm) 

where the symbols q\ and q 2 are the unit vectors parallel to the transverse momenta q 1 
and q 2 , 

and S^(X) is the following factor 

(1) 2-Mn + A) 

n K ' r(n)r(i + A) 

After the integration over the momentum q 2 in Eq.(16) we obtain the product of the 



same vectors [27], |29| 



q ( r ,..^) q ( r ,.,, n ) r(A)r(n + 2A) (2) 

q\ n 2T(2A)r(n + A) n 1 ' 

Then in the front of the expressions, obtained from the contributions singular at e — > 0, 
the following coefficient 

s n (\)=s<»(\)sg\\)- r(n + 2A) 



r(n)r(l + 2A) 



appears. Because S n (\) is factorized and S n (0) = 1, it does not have any influence on our 
results. Thus, the contributions singular at e — > can be calculated easily with the use 



of the formulae from 127, 29 



The ^-independent terms are obtained in a similar way, because one can use of the 
well known formulae (see P9"| ) to calculate the contributions of the Feynman diagrams 
containing ^-functions 

J 1, if 2/ > 
%) -\0, ify<0 

in the integrands. These ^-functions come from an expansion of the expressions 

(li + <?2)~ 2 , (<?i —xq2)~ 2 and ln(gx —q2) 2 in the r.h.s. of eq.(19) (see eq. (|2~4l) , for example). 

However, it is more convenient to work directly with the expansion of propagators 
over the Chebyshev polynomials and to take initially integrals over the angle in the r.h.s. 
of eq.(19). The useful formulae for the expansion of propagators and the properties of the 
Chebyshev polynomials are given below 



ql ~ (& 
(li - Q2) 2 

\n(qi - q 2 ) 2 



00 /a 2 ^ n/2 ' 
l + 2£T„(9i&)(% 

n=l K l2 



\nq 2 2 -2f:-T n (q 1 q 2 )(% 

V?2 



n=l 



n 



0(4 - ql) 

2, «/2n 



qi <-> q2 



d{ql - ql) + 



qi <-> q2 



« - ql 



2T n (x)T m (x) = T n+m (x) + T n _ m (x) (n > m) , 
T n (l) = 1, T (x) = 1, T n (-x) = (-l)"T n (x), 
T 2n (0) = (-1)", T 2n+1 (0) = 0. 



(24) 



(25) 



2.2 The results of calculations 

Applying formulae of the previous subsection to eqs.(16)-(19), we obtain the following 
results for eigenvalues (23): 



6(71,7) 



+ 



n 

7 + - 

1 2 





2n f \ 







2*(l)-*( 7 + ~)-*(l 

X 2 (n,7)-^( 7 + ^)+^(l-7 + |)) 

(?-2C(2)-^)x(».7)-«C(S) 

7r 2 cos (7r 7 ) f/ / nf_\ 2 + 3 7 (1 - 7 ) \ 
1-2 7 )\V V N*J (3-2 7 )(l + 2 7 )y n 



sin (7r 7 )(] 

1 + ^) 7(1 ~ 7) .s 2 

+ iV3;2(3-2 7 )(l + 2 7 ) ^ 



7 + -) -*"(l_ T +-J +2$(n, 7 ) + 2$(n,l- 7 ) 



(26) 



(27) 
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where 5™ is the Kroneker symbol, and *&(z), ty'(z) and are the Euler \1/ -function 

and its derivatives. The function $(72,7) is given below 



$(n, 7 ) 



dx 



x 



rf-X+n/2 



1 + X 



1 v n+l 



- C(2)) + Li 2 (-x) + Li 2 (x) 



00 ^ Xfe 



+ ln(x) ^(n + 1) - (1) + ln(l + x) + £ 

V fc=1 fc + n 

00 fc 



^k + j + n/2 



V'(k + n + 1) - #'(Jfe + 1) + (-l) fc (/3'(A; + n + 1) + /?'(£; + 1' 

+ n + 1) - + r 



(28) 



and 



+ 1 



For the case n = the results (|27|) and ( f2~8| ) coincide with ones obtained in JF). 
Almost all terms in the right hand side of eq. (|27|) except the contribution 



A(n, 7 ) = ( U 



V3 3AU2 



are symmetric under the transformation 7 <-> 1 — 7. Moreover, as it was mentioned in , 
it is possible to cancel A(n, 7) if one would redefine the function g 7_1 by including in it 
the logarithmic factor 



Note also, that the term 
'67 



2C(2) 
9 sv ; 9 JVJ 



10 nA 



x(n,7) = «x(n,7) 



(29) 



is proportional to the LO contribution. It appears often in the MS'-scheme for the pro- 
cesses involving soft gluons and it is absent in the gluon-bremsstrahlung (GB) scheme 
|30|j being more natural for these processes. In the GB-scheme the term ( p9|) can be 
incorporated to a new coupling constant a s (q 2 ) as 



a s (q 2 ) = a s (gV^°)) (ft = yiV c - 1 n f ) 
and eq.(23) is replaced by 



(30) 



& s {q 2 )N c 



7T 



a s (q 2 )N c 

47T 



(31) 
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where 



£(n,7) 



- 6C(3) + 



(|_^)i( xV , 7) _,. (7+i)+ ^ (l _ 7+i 

7r 2 cos(7r 7 ) f / / w/\ 2 + 37(1-7) \ 
in 2 (7r 7 )(l-2 7 )\l + l + N*J (3-2 7 )(l + 2 7 ); 



_ /, , \ 7(1-7) 

I i\r c V2(3-2 T )(l + 2 7 ) " n 



*"( 7 + |) - - 7 + f ) + 2$(n, 7) + 2$(n, 1 - 7) 



(32) 



Analogously to ref. @ we can obtain the eigenvalues (23) in the case of a non-symmetric 
choice of the normalization so of the energy in eq.(15). For example, for the scale sq = q 2 , 
which is natural for deep-inelastic scattering process, we have 



in MS'-scheme 

a s (q 2 )N c 



7T 



X(n,7) + (£(71,7) - 2x{n, 7 )x'(n, 7) 



a a (g 2 )iV c 

47T 



(33) 



and in GB-scheme 

« s (g 2 )iV c 



7T 



X(n,7) + ( 5(n,i) - 2x(n, 7)x'(n, 7) 



«s(g 2 )iVc 

47T 



where 

n 

Considering the limit 7 — > we obtain for n = 



7) = -pxfa, 7) = -*'(7 + ^) + - 7 + £ 



X(0,7) = - + 0( 7 2 ), 5(0, 7 )-2 X (0,7M0,7) = 4 + ~ + C + 

ry* ry 

1(0,7) - 2 X (0, 7 )x'(0,7) = 4 + - + ^ + °(7 2 ), (34) 



where 



.4 
C 



11 , 2n/, 
3 3N* 1, 



n. 



9iV„ 3 



B 



/395 



2C(3) 



yC(2) + ^ 



10JV; + 13j 

(71-18C(2))), 



67 13 n f ' 



(35) 



Following to ref. [0 and using eqs.(31) and (§5|), one can obtain the expression for 
anomalous dimensions of twist- 2 operators 7 at to — + (i.e. near j = 1) 
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in MS-scheme (coinciding with 0) 

a s (q 2 )N c 



7 



7T 

,,3/^2 



+ 



471-3 



1 A 

- + T 

C 



0(w) 



+ 



4n 2 



B 



+ 0(1) 



and in GB-scheme 



7 



a s {q 2 )N c 



7T 



cj 4 



+ 



4tt 2 



+ 



a 3 s (q 2 )N? r 



47T 3 



(36) 



(37) 



Note, that the calculation of terms ~ for anomalous dimensions of the twist-2 
operators at arbitrary values of uj will be performed in the near future (see pi| and 
references therein). 

For the case n > 2 from eq.(|3"3"D one can find also the anomalous dimensions for the 
local operators of the twist t > 2. 



2.3 The asymptotics of cross-sections at s 



oo 



As an example we consider the cross-sections for the inclusive production of two pairs 
of charged particles with mass m in the 77 collision (see ||): 



2 2 1 32/ 



where a em is electromagnetic coupling constant and cr 2 (s) and cr (s) are the coefficients 
for angular dependent and independent contributions, respectively. 

The asymptotic behaviour of the cross-sections Ufe(s) (k = 0, 2) at s — > 00 corresponds 
(for the fixed coupling constant) to the unmoving singularities of /^(t) of the type (a> — 
u;fc) 1//2 , where 



7T 



a s N c 



71 Z 



9tt 5 / 2 



32 J7C(3) (ln(£)) 

7T 5 / 2 



1/2 



^2 



9 • 32^7C(3) - 8 (l n (A; 



l + 0(a SJ j, 

l+0(a s ) 



1/2 



(38) 
(39) 



Here the symbol 0(a s ) denotes unknown a s -corrections to the formfactors. Note, that 
the running of the QCD coupling constant leads to the substitution of the unmoving cut 



by an infinite set of the Regge poles JT2 
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Using our results (26), (27), we obtain the following values for x(k, §) and c(k, \) 
(A; = 0,2): 

X (2 ,h = 4(ln2-l), (40) 



x(o, 

c(0, 



4 In 2, 
11 



2n f \, 67 A/ , 10 n f 1 



22C(3)-64Ls 3 (|) 



If 2 71 - 33^VC(2) 



16 



N 3 

c 



25.8388 + 0.1869^f + 3.8442^ 

N c N* 



(41) 



c(2, 



67 n>/ „ N 10 n f 

— + 2C(2) + f 

9 S V ; 9 N c 



+ 



V3 3iV c y V 

22C(3) + 64Ls 3 (f ) + ^(893 - 3^)^(2) - 64 In 2 



4(ln2 



-3.2652 + 1.5203-^- + 0.3947 

1 V r 



2/ 32 



iVf 



(42) 



where (see [[32| and Appendix A) 

Ls 3 (a;) = -^ln 2 |2sin(|) 

Note, that the function Lss(a;) appears also in calculations of some massive diagrams (see, 
for example, the recent articles |HJ and references therein). 

The LO results (33) coincide with ones obtained in [[|. The value of c(0, 1/2) equals 
to c(l/2) from (see Appendix A). 

As it was discussed in f?J, the NLO correction c(0, 1/2) is large and leads to a quite 
strong reduction of the value of the Pomeron intercept (see recent analyses [|| |TT] of various 
effective resummations of the large NLO terms). Contrary to c(0, 1/2), the correction 
c(2, 1/2) is very small and does not change the small LO value || of the angle-dependent 
contribution. 

In the GB- scheme we have analogously for k — 0, 2 (cf. (31)) 



a s N n 



7T 



—X(k,h 

71 Z 



where 



1 



10 n, 



c(k,-) = c(k,-) + 2C(2) 

Note, that in the case of the non-symmetric choice of the scale so = q 2 in r.h.s. 
of eq.(15), the values of the corrections c(0, 1/2) and c(2, 1/2) are not changed because 
x'(k, 1/2) = 0. Thus, the asymptotics (f41"|) and ([42]) can be applied directly to the deep- 
inelastic scattering process at small values of Bjorken variable x, where x ~ Q 2 /s. 



3 NLO corrections to BFKL kernel in supersymmet- 
ric field theories 



n 



3.1 Gluino and scalar particle pair production 

As it was shown above, the eigenvalue of the BFKL kernel in the next-to- leading approx- 
imation for QCD contains non-analytic contributions proportional to 5° and 5^, which in 
particular leads to a violation of the holomorphic separability. It is natural to investigate 
the BFKL equation in supersymmetric gauge theories where there could be significant 
simplifications. In the case of the N — 1 supersymmetric Yang-Mills theory we should 
consider an additional contribution from a gluino loop. It can be obtained from the QCD 
result for the next-to-leading kernel by a modification of the terms originated from the 
quark loop contributions proportional to n/. Such terms in 5(71,7) contain the factors 
nf/N c oT fif/N^. The terms ~ nf/N c are related to the quark contribution in the running 
coupling constant and will be considered below. 

We discuss here only the terms in 5(n, 7) proportional to rif/N^. Their contribution 
to uj depends on 7 = 1/2 + iv and can be written as follows 



Un f /N§ =-Jn c {-^) U QED , (43) 



where according to eq. (Bl) from Appendix B 

« 2 cos(7r 7 ) / (2 + 37(l-7))ff _ 7(l-7)ffn| \ ( , 

QED sin 2 (7r 7 )(l - 2 7 ) {(3- 2 7 ) (1 + 2 7 ) 2(3 - 2 7 )(1 + 2 7 ) J { ' 

coincides with the position of the Pomeranchuck singularity of the t-channel partial wave 



f w (t) in QED (see |34j]). To give an physical interpretation of this result one should take 
into account, that in the u-plane the BFKL pomeron is the Mandelstam cut originated 
from an exchange of two reggeized gluons. It means, that the amplitude for the transition 
of two reggeized gluons into two others in the t-channel through the fermion loop should 
contain the third spectral function p(s,u). Therefore after the substitution a — > a s it 
should coincide with the amplitude for the virtual photon-photon scattering in QED up 
to the colour factor 



n f ^ tr (t a t b t a t b ) , [t a ,t b ] =if abc t c 



for the corresponding non-planar loop diagram. 

Taking into account, that gluino is a Mayorana particle contrary to the charged elec- 
tron and that its field belongs to the adjoint representation (T a ) bc = —if abc of the gauge 
group, we obtain for the gluino contribution in the iV-extended supersymmetric gauge 
theory the following result 

_ NNj (a s \ 2 

W gluino — TT~^r ^QED ■ (40j 

II \ a / 

Providing that N > 1, there is also a contribution of scalar super-partners of the gluon: 

n s (N) N* (a s y 

^scalar ~ ~ ~ ^SED , (4t>J 

I I \ a J 
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where n s (N) is the total number of scalar particles and ujsed is the position of the 
Pomeranchuck singularity in the scalar electrodynamics (see eq.(B2) from Appendix B 
andref. @): 



« 2 cos(7r 7 ) /(1+tHM , 7(1-7)^| \ 
SED sin 2 (7T7)(l-2 7 ) 1,2(3 -2 7 )(l + 2 7 ) + 4(3- 2 7 ) (1 + 2 7 )J ' 1 J 

For the number of scalar particles in the iV-extended supersymmetric theories we have in 
the case of iV = 2 and N = A 

n,(2) = 2, n a (4) = 6, (48) 

taking into account, that for iV = 4 the gluons and gluinos with both helicities A = ±s 
enter in one super-multiplet with the scalar particles. In particular, the total contribution 
to uj from super-partners of the gluon for iV = 4 is 

Au sp = N 2 [uqed + ^ used) = 

N 2 a 2 cos(7r 7 ) / (11 + 15 7 (1 - 7 )) gj _ 7(1 - 7) ff B | \ 
sin 2 (vr 7 )(l -2 7 ) ^ 4(3 - 2 7 )(1 + 2 7 ) 8(3 - 2 7 )(1 + 2 7 ) J ' 1 J 

It is remarkable, that this term cancels the corresponding gluon contribution non-analytic 
in n (see (23) and 



9 4sin 2 (7r 7 )(l-2 7 ) (3-2 7 )(l + 2 7 )y/ n 2(3 - 2 7 )(1 + 2 7 ) |n| ^ 

3.2 Scalar particle contribution to the gluon Regge trajectory 

Let us consider the contribution of scalar particles to the gluon Regge trajectory in the 
two- loop approximation. We use the results of ref. |2"2" ], where the corresponding quark 
correction was calculated, which gives a possibility to find easily the gluino correction 
using the substitution rif — > N N c . The quark contribution can be presented in the form 



n 9 (g) 



2(1 + e) r 2 (i + e)r(i-£) 

qKH ' ~ (3 + 2£)(l + 2 £ ) £ r(l + 2£) (q>)- £ ' 

where all (j (g) (g^q 2 — q^q^) is one-loop quark self-energy correction to the polarization 
tensor, and a the is bare QCD constant 

_ g 2 N c T(l - e) 
a (4vr) 2 + £ 

with well-known relation between g 2 and g 2 in MS scheme: 



g 2 = g>~ 2£ 



1 /ll 2nt . 

1 + i(y-3ivf^ 



(51) 
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where the last term in the r.h.s. of Eq.flSTD is proportional to the leading order term f3$ 
(0) of QCD ^-function. 

The scalar particle correction can be presented in the similar form 

«2>(*) = "^P • ^ / [^sM ~ IW,)] , (52) 

n sc (g) = 



3N C 2 J q 2 (q-q) 2 

1 r 2 (l+e)r(l-e) 
"2(3 + 2e)(l + 2e)er(l + 2e) ( q 2)~ e ' 



where all sc (g) (g^q 2 — q^q u ) is one-loop self-energy correction to the polarization tensor 
from scalar particles and a is the bare constant of N = 4 supersymmetric theory 



_ 9 2 N C 
a 



(An) 2 

As it is well known, the coupling constant of the N = 4 supersymmetric theory does not 
run, i.e. its /3-function is zero, and, hence, a = a^. 



3.3 The results of calculations 

Putting all things together, in the case of N = 4 supersymmetric theory we have the 
following representation for the Regge trajectory 



susy 



(q) = -a,. 



a l 



2C(2) 



+ 2 In 
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2C(3) 



(53) 



Using (53), by analogy with the results of [0] for QCD we can easily find the corre- 
sponding kernel K^ usy (qi, q 2 ) for the N = 4 SUSY 



Kr v (qi,Q2) 



-2c 



7r 1+£ r(l - e) { ( qi - q 2 ) 



+a. 



' (qi - 92) 



_-2C(2) + e^-- + 14C(3) 



In 



+a. 



(qi - <?2) 2 (<?i + q2) 2 

{ql - ql) 



{ql + ql) 4 



-(-§: 



-a: 



1 



(?i -q2) 2 (qi + q2) 2 , 



dx 



1 J (q 2 - xqi) 2 



In 



' xq\ 

A. 



(54) 



The set of eigenvalues for eigenfunctions (22) of the homogeneous BFKL equation in 
the framework of N = 4 supersymmetric theory 



susy 



4a 



X(nn)+S msy (n n )a 



(55) 
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can be found analogously to the previous section in the following form 



S susy (n,j) 



2$(n, 7) + 2$(n, 1 - 7 ) - ( i - 2((2) ) x (n, 7) 



6C(3)-*"(7 + f)-*"(l-7 + f) 



(56) 



where the function $(n, 7) was defined in eq.(28). 
Note, that the term 

1 -2C(2)) X (n,j) 



appears as a result of the use of the MS-scheme in intermediate calculations (see (|29|)). 
This term can be eliminated by a redefinition of the coupling constant 

1 



a->a = a+l^- 2C(2)Ja 2 , 

which is equivalent to the use of the BG-scheme (see (|30~D ). For the new coupling constant 
a above results can be written in the following form 



to 



susy 



4a 



x (n,7) + ^(n,7)a 



where 

~5 susy (n^) 



2$(n, 7) + 2$(n, 1 - 7) - 6C(3) - *"( 7 + ~) - 



n 

7 + 7T 



• (57) 



Analogously to the previous section and ref. [|7| we can obtain the eigenvalues of the 
kernel in the case of non-symmetric choice of the energy-normalization sq in eq.(15). For 
the scale s = q 2 we have 



in MS-scheme 



UJ 



susy 



4a 



X (n, 7) + 5 sus «(n, 7) - 2 X (n, 7 )x'(n, 7) a 



(5f 



and in GB-scheme 



a; 



susy 



4a 



X(n, 7) + o" suss/ (n, 7) - 2 X (n, 7 )x'(n, 7) 5 



The anomalous dimensions 7 susy of the twist-2 operators at uj — > (and n = 0) are 

- J^susy 



in MS'-scheme 



y«n* = 4 _ 

and in GB-scheme 

7 sns ^ = 4 5 



- + 0(u] 

UJ 



UJ 



' /~<susy 

+ 0(l))+a 2 (^ + 0(^- ] 



UJ 



(_ + 0(w) ) +s (_ + 0(1) ) +s .( 



(Jsusy 



UJ 



+ uj 



-1 



(59) 



(60) 
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where 



B susy = -, B susy = 2C(2), and C susy = 2((3). (61) 
3 

For the case n > 2 one can calculate also the anomalous dimensions of the correspond- 
ing local operators. Because in the case N = 4 SUSY the result is analytic in |n|, one can 
continue these anomalous dimensions to the negative values of \n\. It gives a possibility to 
find the singular contributions of the anomalous dimensions of the twist- 2 operators not 
only at j — 1 but also at other integer points j = 0, —1, —2.... In particular, in the Born 
approximation for the anomalous dimension of the operator corresponding to the sum of 
distributions of gluons, gluino and scalar particles we obtain 7 = 4a(\l/(l) — ^(j — 1)) 



which coincides with the result of the direct calculations (see [^0j). Thus, in the case 
N = 4 the BFKL equation presumably contains the information sufficient for restoring 
the kernel of the DGLAP equation. 



4 Conclusion 

In this paper we calculated next-to-leading corrections to eigenvalues of the BFKL kernel 
in QCD for arbitrary conformal spins n (see (p7[)). Now with the use of eqs. (39), 
(42) one can find the spin-dependent contribution to the total cross-section for the 7*7* 
annihilation into hadrons at high energies. We calculated also next-to-leading corrections 
to the BFKL equation in the supersymmetric field theories and next-to-leading corrections 
to anomalous dimensions of twist-2 operators (see (56), (59)). In a particular case of the 
extended N = 4 SUSY the result of our calculations is significantly simplified (see eqs. 
([36]) and (|57|)). The absence of the coupling constant renormalization in this model 
leads presumable to the Mobius invariance of the BFKL equation in higher orders of the 
perturbation theory. This invariance can allow one to guess the next-to-leading corrections 
to the BFKL kernel in QCD for a general case t 7^ 0, which is very important for the theory 
of the processes at small x. The remarkable cancellation of non-analytic contributions 
proportional to and 5% at N = 4 SUSY is a possible manifestation of the integrability 
of the reggeon dynamics in the Maldacena model |p2TJ corresponding to the limit iV c — > 



00. Note, that in this model the eigenvalues of the LLA pair kernels in the evolution 
equations for the matrix elements of the quasi-partonic operators are proportional to 
— 1) — pOfl , which means, that the corresponding Hamiltonian coincides with 



the local Hamiltonian for the integrable Heisenberg spin model. The residues of these 
eigenvalues in the points j = —k can be obtained from the BFKL equation by the analytic 
continuation of the anomalous dimensions to negative integer values of the conformal spin 
Therefore the DGLAP equation is not independent from the BFKL equation in N = 4 
SUSY and their integrability properties at N c — > 00 are related. We shall return to these 
problems in our future publications. 
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A Appendix 



Here we demonstrate how the function Lss(x) appears in eqs. (|41|) and fl42] 
The direct evaluation of the values of c(0, 1/2) and c(2, 1/2) leads to 



c(0,i) 



c(2,i) 



,11 2n f V 67 > , , 10 n f 1 

2 *- ln2 + 2C 2 + + 

1 3 3 NJ 9 su 9iV c 41n2 



22C(3) 



+ 16 / — arctan(\/t) In 
Jo t 



1 \ l_ 
1 -t) + 16 



■S] 1 33^VC(2) 



(A.l) 



,11 2n f \. s 67 A/ , 10 n/ 1 



22C(3) 



- 16 ^ 1 ^axctan(Vt)hi( r ^)+^(l89-3^)7rC(2)-641n2 



Note that ( |A.1|) coincides exactly with eq.(16) of 0. 

To obtain eqs. ( Ell ) and (K2J) we should evaluate the integral 

/ — arctanfv t) m 

Using the following substitutions 

2 / ■ \ i . 1 + y — ix 

t = p , arctanfx + ly) = - In h Ink 

2 1 — y + ix 

we obtain for it the following representation (for k = 0) 



3 i7T/2>, 



(A.2) 



/ — arctan(v / t) In ( —^—) = -i [ — fln(l - p) + ln(l + p) 
io t VI — i/ Jo p \ 



where [55 1 



- t 

ln(l — ip) — ln(l + zp) 



F 1 (i) + F 2 (i)-F 1 H)-F 2 (-i) 



(A.3) 



Fi(a) = / — ln(l — p) ln(l + ap) = Li 3 (— a) + Si 2 (— a) and 
jo p 

^ ln(l + p) ln(l + ap) = S M - S 1>2 + S M (-a) 



F 2 fa) 



o p 
+ Li 3 (a) 



-C(3) + -ln 2 aln 



1 + a 
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— In a 



C(2) - Li 2 



1 — a 
1 + a 



(A.4) 



Here S ntP (x) are the Nielsen generalized polylogarithms (see [3E 



/ «y ln n-i 2/ i n (i_ X2/ ) 



(n — l)!p! jo y 



(A.5) 



and the polylogarithms Li n (x) = S n _i i i(z). 
With the use of (|A.4f) we have 



Fx(±i) 
F 2 (±i) 



Li 3 (^i) + Si, 2 (=r), 

2S 1>2 (^) - S 1)2 (±i) + Lia(±i) - ^C(3) ± 

o Z 



Ii 2 (=Fi) + ^C(2) 



(A.6) 



Using the integral representation (|A.5| ), we have for the Euler dilogarithm Li 2 (i): 

dp 



Li 2 (l) - Li 2 (z) = 
If we replace s by s = e tB corresponding to 



l- s =l- e 



2 sin ~ e 1 ^' 2 



we obtain 



where 



Li 2 (l)-Li 2 (z) = ^C(2)-iG, 
In a; 



G = - / dx 



1 + X 1 



Cl 2 (vr/2). 



Here according to ref. [32 



I In"" 1 


2sin(|) 


/o 





dy. 



Thus, 

Li 2 (±z) = -^C(2)±iG. 
Repeating above calculations, one can derive the relations 
Li 3 (±i) 



-— C(3)±— «7rC(2), 
32 sv ; 16 SV ; ' 



Si )2 (±i) 



6l C(3) "4 GT 2 



Ls 3 (0 + ^C(2) 



The obtained results flPD-flOD are in agreement with ref. |32l 
Putting QA77D-(pO|) to (pOQ and QO|) , we have 



r 1 dt / l 

/ — arctanfv £) In 

Jot V 1 — t 



4Ls 3 ( - J - yvrC(2) 



and, thus, eqs.( |A.l|) and ( |A.2j ) coincide with Eqs.(£|l]) and (|42T) , respectively. 
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(A.7) 

(A.8) 
(A.9) 



B Appendix 

Although the Pomeranchuck singularity in QED was investigated many years ago (see 
||34j|), here we review the obtained results and generalize them to the case of the charged 
scalar particle production. 

To begin with, let us consider the electron-electron scattering in the Regge kinematics 
%PaPb = s ^> q 2 , where pa and ps are momenta of the colliding electrons. In the Born 
approximation the corresponding amplitude equals 

5^ u 2/s s 

A ee ^ ee = j^p A ) —j v (p B ) = j (p A )p B —j v (p B )p v = 2e 2 — , 
q z q z q z 

where j^ip) = e u(p') 7^ u{p) . The elastic cross-section is 

1 f d 2 q \A ee ^ ee \ 2 2 f d 2 q 
a ei = - / 77; Tn" 5 = 4 « 



4 7 (2tt) 2 s 2 J |g| 4 ' 

This cross-section is divergent at small q due to an infinite radius of the Coulomb in- 
teraction. We shall ignore this circumstance because in the case of the scattering of 
non-charged particles A and B the integrand is multiplied by a product of their impact 
factors 5>a(<?) and §B(q) vanishing at q = 

d 2 q 



VAB = J—r^M<l) < S>B(<l) 



(2n) 2 \q 

which leads to the cancellation of this divergency (see pl| ). Note, that for the electron 
being a charged particle the impact factor is non-zero for q = 0: 

$ e (g) = Area 

The amplitude for the production of a massless electron-positron pair in the ee scat- 
tering is 

A - 9p 4 _!_ A 

JT-ee — >eee~ t ~e 2 2 7*7* — >e + e 5 

<?1<?2 

where the factor A~* 7 *^ e + e describes the transition of two virtual photons into the e + e~ 
pair 

A^^e+e- = u e -(h) A e v e +(-k 2 ) , qi - q 2 = h + k 2 . 
The amplitude A e is given below 



where we introduced the Sudakov variables for momenta k\ and k 2 of the produced electron 
and positron 
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k r = a r pB + PrPA + kr , sa r (3 r = —{k^) 2 = \k r \ 2 
and used for the momenta of the virtual photons the following relations 

qi = (Pi + AO Pa + qi , q-2 = -(ai + a 2 ) Pb + q 2 

valid at high energies when 2k\k 2 <C s (corresponding to small a r and f3 r ). Note, that 
with the use of the Dirac equation for the produced particles we can rewrite A e in the 
other form 

A _^ 2 / sfcqj + qj qi Pb _ s_Pi gj + Pb q 2 Qi 
e ~* s{(3 1 + (3 2 ) [ik.-q.f+lhl 2 ^ \k 1 + q 2 \ 2 + \k 1 -q 1 + q 2 \ 2 ^ / 

and verify its vanishing at q± = or q^ = 0. 

Analogously for the production of a pair of charged scalar particles s we have 

A , - ?r 4 S A 

mi 

where 

A s = 2 + — -7} — ,9 „ h 



2 


N' 


2 


\h - gi| 


2 + 


\h\ 


2 §2 
Pi 



2 


^2 


2 


&2 - qi\ 


2 + 




2 & 

ft 



It can be written in the form 

2 / 2& - g?) Pi(ql,2k 2 -qt) 



A 



S P1+P2 \\k 1 -q 1 \ 2 + \k 1 \ 2 ^\k 2 -q 1 \ 2 + \k 2 \ 2 % i 

to make obvious its vanishing at q\ = or at q 2 = 0. The total cross-section for the pair 
production in the collision of the particles A and B after calculating the integrals over ctj 
(i — 1, 2) with the use of the (^-functions 5(saii(3i — \kif) can be presented as follows 

(Tpair = log -J / 4 7^ ^ 2 | |4 ^1 > ?2 ) • 

m J (2ti) 2 \qi\ {2n) 2 \q 2 \ 
The kernel K(qt,q^) for electrons e and scalar particles s is correspondingly 

K e , s {qx, 82) = 16a 2 / / rfa; / e ,,(fct; x) , x = ^ 

J {2ir) z Jo pi + p 2 

and 



P tr{k x A e k 2 At) f(V \A S 
fe(h;x)= — — , f s (ki;x) = 



16x(l — x) ' ' lQx(l — x) 

The factor log(s/m 2 ) was obtained as a result of the integration over /3 — f3\ + f3 2 in 
the region m 2 / s <C f3 <C 1, where the mass scale m equals to an essential value for 
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transverse momenta of charged particles inside the colliding particles. We obtain for the 
e + e~ production 

t (t* \ \\h\ 2 \h - q\\ 2 \\k 2 \ 2 \h + q 2 \ 2 

fe{ki;x) = - * — -2 + 



(jki - xq^ 2 + x(l - x) |gi| 2 ) (\h + q 2 - xq^ 2 + x(l - x) |g 2 | 2 ) 

(fei - qt, q2 + fei) (fa, fei) - (fci - gi , fa) (fci + ?2 , fci) - (fci + fg, fcgj (fci - gi , fci! 
(|fci -xgi| 2 + -x) |gi| 2 ) + g 2 -xgi| 2 + - x) |gi| 2 ) 

and for charged scalar pair production 

, f-r* x \x(l-x)(qt,2k[-qt) 2 \ x{l - x) {qt, 2k* 2 - qt) 2 

fs(k 1 ;x) = - r i -2 + 



(jki - xqi\ 2 + x(l - x) |gi| 2 ) (\h + q 2 - £<?i| 2 + x(l - x) \q 2 [ 

| x(l - x) {qt, 2kl - qt) (qt, 2%, - qt) 

([fa - xq^ 2 + x(l - x) |gi| 2 ) (\h + q 2 - xq x f + x(l - x) |gi| 2 ) 

By introducing the Feynman parameter y to unify two factors in the denominators of 
the last contributions to f e (k\;x) and / s (k\,x) and integrating the kernel over fci one 
can present it in the form (cf. ||34j| ) 

_> 8a 2 f 1 f 1 _> _> 

K e ,s{qi,q 2 ) = / dx dy ip etS (q 1 ,q 2 ;x,y) , 

7T JO JO 

where for the electrons 

<Pe(qi, Q2; x, y) = 

(x(l -x)+ x(l -y)- 4a; (1 - x)y(l - y)) {q^ 2 \q 2 \ 2 - 2x(l - x)y(l - y) (qt, q^f 

x(l -x) \qi\ 2 + y(l-y) \q 2 \ 2 

and for scalar particles 

(x(l -x)+ x(l -y)- 8x(l - x)y(l - y)) |gi| 2 \q 2 \ 2 + 4x(l - x)y(l - y) (qt, g^f 

A(x(l-x)\ qi \ 2 + y(l-y)\q 2 \ 2 ) 

Note, that for the supersymmetric QED, where there are scalar charged particles, the 
total contribution to the kernel is 

<Pe(qi, Q2; x, y) + 2(p s (qt, <gj x, y) = 3x(l - x) + 3x(l - y) - 16x(l - x)y(l - y) 
\Qi\ 2 \Q2\ 2 2(x(l-x)\ qi \ 2 + y(l-y)\q 2 \ 2 ) 

and does not contain the term proportional to (qt , q^) 2 '■ 

The total cross-section in LLA can be written in the form 

ra+ioo duo { s y f d 2 q 1 $a(<?i) d 2 q x <$> B (q 2 ) , 
atot -l ■ 2^U*J J (2^\ qi \ 2 (2n)i\q 2 \ 2 



a— loo 
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where the t-channel partial wave f^(oi,Q2) for the 7*7*-scattering satisfies the Bethe- 
Salpeter equation 

= (^) 2 5 2 (q -q')+ [ .f. q \ 2 2 • 
J (2tt) 2 

The solution of this equation can be obtained in the form of the Fourier expansion 

o /I \\2iv 1 

2,77 ( \Qi\\ 1 



_J-oo fli 02 V (79. 



9i| |92 1 V | <72 1 y w - w(n, i/) ' 



where is the phase of the two-dimensional vector ~ql written in the complex form as 
\qk \ e l<Pk and u(n, v) is the eigenvalue of the homogeneous Bethe-Salpeter equation: 

•/ (27r) 2 |g| Vkly 

Using the above expressions for K e s (lf, q), it is easy to calculate u(n,u) for the 
spinor and scalar electrodynamics (cf. 



a 2 sinh^) 1 1 t 1o ^ x0 l+4z/ 2 

32 (cosh(7rz/)) 2 1 + v 2 v 



w Q En(n,v) = - 11/1 ^ , -J- ± (ii + 12 , 2 ) - ft, , (B.1) 



The high energy behaviour of the total cross-section is 

oo 

Otot ~ E c ™ sA(n) c °s(n(v?' - V 2 )) 

ra=0 

where A is the Pomeron intercept: 

A QED (n) = rra 2 £ _ J_ (B.3) 



A 5BD (n) = vra 2 5° n + — 5 2 n j . (B.4) 

Note, that for the supersymmetric QED the results are especially simple: 

. . a 2 sinh(7rz/) rn , . 7ra 2 rn 

ususrin, f ) = — - — — — rra 6 n , AsusY{n) = —8 n . 

£V (cOSh(7Tl/)) I 
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